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D . Abstract 

Q 

' We show that Einstein equations are compatible with the presence of massive 

' point particles and find corresponding two parameter family of solutions. They are 

defined by the bare mechanical mass Mq > and the Keplerian mass M < Mq of 
I the point source of gravity. The global analytical properties of these solutions in 

QQ ■ the complex plane define a unique preferable radial variable of the problem. 

00 ■ 

O 

^ ■ 1 Introduction 

^ I A clear physical motivation for consideration of massive point particle sources of gravita- 

tional field in GR, both electrically neutral and charged ones, can be found as early as in 
1962-63 Feynman lectures on gravity jl]. In spite of this fact the problem is still open, 
bi). At present the vast majority of relativists do not accept the consideration of point par- 

ticles in general relativity (GR) as incompatible with Einstein equations (EE) idealization. 
There are different reasons: some doubts about consistence of the theory of mathematical 
' distributions (like Dirac (5-function S{r)) with the nonlinear character of EE; the clear un- 

■ derstanding of the drastic change of geometry of the Riemannian space-time Wl^^'^^{g^iy} 

in a vicinity of a point with infinite concentration of energy in it, etc. 

On the other hand it is obvious that in Nature very distant objects like stars look 
like "points" of finite mass and finite luminosity. This fact has a proper mathematical 
description in Newton theory of gravity in the language of mathematical distributions. 
A formal mathematical problem is to find a corresponding idealized treatment of such 
objects in GR, as well, but up to recently no reasonable approach was known. 

Here we show that correct mathematical solutions of EE with 6{r) term in the rhs do 
exist. Such solutions describe a two parameter family of analytical space- times M(i'3){^^^} 
with specific strong singularity at the place of the massive point source with bare me- 
chanical mass Mq > and Keplerian mass M < Mq. 

The price, one has to pay for this enlargement of standard GR framework, is: 
1) To consider metric coefficients g^i/{x) as functions of class C*^ of the coordinates x. 
Some of them are to have definite finite jumps in their first derivatives at the place of 
the point source, necessary to reproduce the 6{r) term in the Ihs of EE from the Einstein 
tensor G"^ ~ (^(i")- 
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2) To accept the unusual geometry of space-time around the matter point in GR 
which actually appeared at first in the original Schwarzschild article |2| and has been 
discussed by Brillouin j3] as early as in 1923. This geometry is essentially different from 
the geometry around space-time points with finite energy density in them. The global 
properties of the space-time are different from the ones of popular Hilbert-Droste-Weyl 
form of the Schwarzschild solution |i4j. 

At present the original Schwarzschild geometry and other similar geometries of space- 
time are widely ignored in GR. In addition, the historical facts unfortunately are not 
reproduced in completely correct way in the most of the modern literature on this subject. 

The important physical consequences, which can be derived using the new solutions 
of EE, discussed in the present article, can be considered after these solutions will be 
properly studied. 

2 Formulation of the Mathematical Problem 

In its proper frame the single point particle with bare mechanical mass Mq, placed 
at the origin of the standard spherical coordinate system in the 3D Riemannian space 
M^^^{gij} C M^^'^^{g^^} yields the familiar static metric [5j 



with three unknown functions gtt{r) > 0, gr-rir) < 0, p(r) of the radial variable r > 0. The 
variable r is not defined by the S0{3) symmetry of the problem, nor by its global-time 
translation invariance with respect to the group Tt{l). The only clear thing is that the 
value r = must correspond to the center of spherical symmetry, C, where the massive 
matter point is placed. 

In contrast, the quantity p has a clear geometrical meaning: it defines the area Ap = 
Airp"^ of the sphere, centered at the point C. Physically p defines the luminosity of the 
objects, which is reciprocal to Ap. Therefore we shall refer to p as "a luminosity variable" . 

At first glance the function p(r) may be chosen in quite arbitrary way, thus fixing the 
remaining (radial) gauge freedom of the problem - the only one, which is not fixed by 
symmetry reasons. In this article we shall show that choosing a definite class of functions 
p(r) one can solve correctly the Einstein equations: 



with stress-energy tensor ~ Mo5(r), which describes a massive point source with bare 
mass Mo and mechanical action Amo = —Mq J ds. 

It may seems strange that to solve the problem, one needs to fix the class of coordinates 
by choosing the radial gauge. As we shall see, the proper choice of coordinates is dictated 
naturally by the boundary conditions. In the problem at hand these conditions are masked 
in Dirac 5-function, which describes in a formal mathematical way the properties of source 
of gravity and its boundary. 

The shortest way to derive the field equation of our specific problem is to perform a 
restriction on the orbits of the group S'0(3) x T((l) of the total action Atot = AGR + Amq, 
where Agr = —j^^q- J d'^x^yJg\R is the Hilbert-Einstein action. Neglecting the inessential 



ds^ = gtt{r) dt^+g„{r) dr'^ - p{rf{d9'^ + sm 
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surface terms and using r as an independent variable, we arrive at one dimensional problem 
with Lagrangian 
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Here Gn is the Newton gravitational constant, 6{r) is the ID Dirac (5-function (We 
are using units c=l.) The corresponding Euler-Lagrange equations read: 



where the symbol " = " denotes a weak equality in the sense of theory of constrained 
dynamical systems. (As a result of the rho-gauge freedom the field variable y/—grr is not 
a true dynamical variable but rather plays the role of a Lagrange multiplier, which is 
needed in a description of constrained dynamics.) Of course, the same equations ()2.2|) 
can be derived directly from EE ()2.H) . 

Outside the point source of the gravitational field (i.e. for r > 0, where 6{r) = 0) one 
obtains the standard solution of this system [7j: 



where pc = 2GnM is a new integration constant - the Schwarzschild radius, M is the 
gravitational (Keplerian) mass of the source, and p{r) is an arbitrary function of r. 

As usual, this solution describes the gravitational field outside the point source (and 
outside any body with spherically symmetric mass distribution). According to the 
Birkhoff theorem jH] it is unique (up to the choice of the radial variable r). Nevertheless, 
being a pure vacuum solution of EE, it is not able to describe the mass distribution Mo5(r) 
in the very point source, as well as the mass distribution in any other matter source. 

To see this in a more transparent way, let us consider the Hilbert gauge (HG) of the 
radial variable: r = p in which the field equations ()2.2|) can be rewritten in the form: 




(2.2) 



gtt{r) = 1 - pa/pir), gtt{r) ^„.(r) 



(2.3) 
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Here we have utilized the properties of Dirac 5-function and the assumption that the 
metric coefficients are usual functions of r (not distributions). Then the only remnant 



of the function p{r) in the system Eq. ()2.4j) are the numbers po = p(0) and ao = sign{p'{0)), 
which enter only the first equation, related to the source of gravity. 

One usually ignores the general case of an arbitrary value Po 7^ accepting the value 
Po = 0, which seems to be natural in HG. Indeed, if we consider the luminosity variable as 
a measure of the real geometrical distance to the point source of gravity in the 3D space 
(which is not the case), we have to accept the value po = for the position of the point 
source. Otherwise the (5-function term in the Eq. ()2.4|) will describe a shell with radius 
Po 7^ 0, instead of a point source. Such conclusion is based on the above interpretation of 
the variable p. 

Actually the point source has to be described using the function 6{r) and is placed 
at the point r = by definition. There is no reason to change this original position of 
the source, or the interpretation of the variables in the problem at hand. To what value 
of the luminosity variable po = p(0) corresponds the real position of the point source is 
not known a priori. This depends strongly on the choice of the rho-gauge. One can not 
exclude such a nonstandard behavior of a physically or mathematically reasonable rho- 
gauge function p(r) which leads to some value po 7^ 0. Physically this means that instead 
to infinity, the luminosity of the point source will go to a finite value, when the distance 
to the source goes to zero (see, for example, 0). This interesting possibihty appears in 
curved space-times due to their unusual geometrical properties and is not supported by 
our Euclidean experience. If one accepts the value po = 0, one has to recognize that 
the HG singularity at p = will be space-like, not time-like, because for p G [0, pc] the 
variables p and t are changing their roles. This will be a quite unusual and non-physical 
property for physical source of a physical field of any kind. 

The solution of the subsystem formed by the last two equations of ()2.4|) is well known: 

gu{p) = l-pG/p, 9pp{p) = -^/gttip) (2.5) 

- the Hilbert form of the Schwarzschild solution. Note that in this subsystem one of the 
equations is a field equation, but the other one is a constraint. However, these functions 
do not solve the first of the Eq. (j2.4p for any value of po, if Mq ^ 0. Indeed, for these 
functions the left hand side of the first field equation equals identically zero and does not 
have a (5-function-type singularity, in contrast to the right hand side. Hence, the first 
field equation remains unsolved by the functions ()2.5p . Thus we see that the assumption 
that gttij), grrij) and p(r) are usual & smooth functions, instead of distributions, yields 
a contradiction, if Mq 7^ 0. Indeed, in this way one is not able to describe correctly the 
gravitational field of a massive point source of gravity in GR. 

This mathematical result is the real basis of the widespread opinion, according to 
which it's impossible to describe a massive point in GR. Actually, the right conclusion is 
that the class of (and even the class of C^) metrics is not proper for this purpose. To 
reach such goal, the first derivative with respect to the variable r of at least one of the 
metric coefficients g^y must have a strictly definite jump. It is needed to reproduce the 
Dirac (5-function in the energy-momentum tensor Tj^ of the massive point particle via the 
Einstein tensor C^. Otherwise the EE ()2.1|) can not be fulfilled. 

The above consideration illustrates one more important juncture. It is obvious that 
physical results of any theory must not depend on the choice of the variables and, in 
particular, these results must be invariant under changes of coordinates. This requirement 
is a basic principle in GR. It is fulfilled for any already fixed mathematical problem. 



Nevertheless, the change of the interpretation of the variables may change the formu- 
lation of the mathematical problem and thus, the physical results, because we are using 
the variables according to their meaning. For example, if we are considering the lumi- 
nosity variable p as a radial variable of the problem, it seems natural to put the point 
source at the point p = 0. In general, we may obtain a physically different model, if 
we are considering another variable r as a radial one. In this case we shall place the 
source at a different geometrical point r = 0, which now seems to be the natural position 
for the center C. The relation between these two geometrical "points" and between the 
corresponding physical models strongly depends on the choice of the function p(r), i.e. 
on the radial gauge. Thus, applying the same physical requirements in different "natu- 
ral" variables, we arrive at different physical theories, because we are solving EE under 
different boundary conditions, coded in corresponding Dirac 5-functions. One has to find 
a theoretical or an experimental reasons to resolve this essential ambiguity. 



3 Field Variables and Radial Gauges, Suitable for a 
Correct Treatment of the Point Source of Gravity 



An obstacle for the description of the gravitational field of a point source at the initial 
stage of development of GR was the absence of an adequate mathematical formalism. 
Even after the development of the correct theory of mathematical distributions [U] there 
still exist an opinion that this theory is inapplicable to GR because of the nonlinear 
character of Einstein equations ()2.H) . see for example |H]. In a recent article the 
authors have considered singular lines and surfaces, using mathematical distributions. 
They have stressed, that "there is apparently no viable treatment of point particles as 
concentrated sources in GR" . 

Here we propose a novel approach to this problem, based on a specific choice of the 
field variables in the metric: 



where V'i(r), f2{r) and (f{r) are unknown functions of the variable r and p is a constant - 
the unit for luminosity variable p = p e"'^^"'"'^^. The corresponding form of the Lagrangian: 



shows that the field variables v^i(r), (p2{r) and ip{r) diagonalize the kinetic-like part in L. 
Hence, they play the role of a normal fields' coordinates in the problem at hand. Now 
the field equations read: 
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(3.3) 



where Ar = e '^-^ (^^^) is related to the radial part of the 3D-Laplacean. The variation 
of the total action with respect to the auxiliary variable gives the constraint: 



(3.4) 



The advantage of the above normal fields' coordinates is that when expressed in them the 
field equations ()3.3p are linear with respect to the derivatives of the unknown functions 
(y9i_2(r). This circumstance legitimates the correct application of the mathematical theory 
of distributions and makes our normal coordinates privileged field variables. 

The choice of the function <f{r) fixes the rho-gauge in the normal coordinates. We have 
to choose this function in a way that makes the first of the equations ()3.3|1 mathematically 
meaningful. Note that this inhomogeneous equation is quasi-linear and has a correct 
mathematical meaning if, and only if, the condition |</'i(0) — v^(0)| < oo is satisfied. 

Let's consider once again the domain r > 0. In this domain the first of the equations 
flH.Hj) gives 93i(r) = Ci J e~'^^^^dr + C2 with arbitrary constants Ci,2- Suppose that the 
function (p{r) has an asymptotic exp(— (^(r)) ~ kr^ in the limit r +0 (with some 
arbitrary constants k and n). Then one easily obtains ipi{r) — i^{r) ~ CiA;r"+^/(n + 1) + 
n In r + In k + C2 - a n ^ —1 , and (pi{r) — ip{r) ~ {Cik — 1) Inr + In A; + C2 - for n = —1. 
Now we see that one can satisfy the condition lim^-^o ~ = constant < 00 for 

arbitrary values of the constants Ci^2 if, and only if n = 0. This means that we must have 
(^(r) ~ A; = const 7^ ±00 for r — >■ 0. We call such gauges regular gauges for the problem 
at hand. Then (pi{0) = const 7^ ±00. Obviously, the simplest choice of a regular gauge is 
(^(r) = 0. Further on we shall use this basic regular gauge (BRG). Other regular gauges 
defer from it by a regular rho-gauge transformation which describes a diffeomorphism of 
the fixed by the BRG gauge manifold 'M^^^ {gmnii^)} ■ In terms of the metric components 
the BRG fixing condition reads p'^gu + p^grr = 0. 

Under this gauge the field equations ()3.3|1 acquire a simple quasi-linear form: 

^lir) = e^^W 5(r), y,^'(r) = 1 e^^^^ (3.5) 

and constraint (|3.4|) is: 

-(pVi)' + (pV2)'-e'^^ = 0. (3.6) 

Hence, the BRG gauge (^(r) = has the unique property to split the system of field 
equations (|3.3|) and the constraint (|3.4|) into three independent relations (j3.5p . (j3.6|) . 



4 Solution of Einstein Equations with Massive Point 
Source 

4.1 Solution in the Basic Regular Gauge 

The new form of Einstein field equations ()3.5p of our problem has a simple general solution 
in terms of mathematical distributions (called sometimes "generalized functions" 0): 

¥^i(r) = ^|^e'^^(°)(e(r)-0(O))r + y.;(O)r + v;i(O), 

<^2(r) = -ln(^-i=sinh(^V2i;^;^2e^^^ . (4.1) 

The first expression in Eq. ()4.1|) represents a distribution (pi{r). In it G(r) is the Heaviside 
step function. Here we use the additional assumption 0(0) := 1. It gives a specific 



regularization of the products, degrees and functions of the distribution 0(r) and makes 
them definite. 

The second expression (p2{r) in Eq. ()4.H) is a usual function of the variable 
r. The symbol r^o is used as an abbreviation for the constant expression r^o = 
sign (v?2(0)) p sinh (v^e-'^^Co) j /^/2e^. The constants 

., = -ipVi(r)V^^;(O)e-(°)(e(r)-0(O)), e, = \ (pM(r)^-e^-«) (4.2) 

are the values of the corresponding first integrals (j4.2j) of the differential equations (|3.5|) 
for a given solution (j4.H) . 

Then for the regular solutions ()4.1|) the condition ()3.6|) reads: 

^i + ^2 + ^|^e^^W(e(r)-e(0)) =0. (4.3) 

An unexpected property of this relation is that it cannot be satisfied for any value of 
the variable r G {—oo, oo), because ei^2 are constants. The constraint ()4.3|1 can be satisfied 
either on the interval r G [0, oo), or on the interval r G {—oo, 0). If, from physical reasons 
we chose it to be valid at only one point r* G [0, oo), this relation will be satisfied on the 
whole interval r G [0, cxd) and this interval will be the physically admissible real domain 
of the radial variable. Thus one can see that our approach gives a unique possibility to 
derive the admissible real domain of the variable r from the dynamical constraint ()3.6p . 
i.e., this dynamical constraint yields a geometrical constraint on the values of the radial 
variable. As a result, in the physical domain the values of the first integrals (j4.2j) are 
related by the standard equation 

Stot = €1+62 = 0, (4.4) 

which reflects the fact that our variation problem is invariant under local re- 
parameterization of the independent variable r. At the end, as a direct consequence 
of relation ()4.4|) one obtains the inequality £2 = —€i > 0, because in the real physical 
domain r G [0, 00) we have ei = —^p^i^\{r)^ = const < 0. 

For the function pBRcij) > 0, which corresponds to the BRG, we obtain 

PBRG{r) =Pg{i- exp (4^^— p^)^ . (4.5) 

Now we can fix the arbitrary integration constants v'i(O), V5'i(0), and €2'- 

1) Imposing a number of additional requirements: i) to have an asymptotically flat 
space-time, taking into account that in the present variable r we have p(r) 00 when 

^ Tco'i ii) to have a correct Keplerian mass M for p — 00; iii) to have a consistence with 
the relation gttQrr + p'^ = 0; 

2) Making use of suitable choice of the units for the luminosity variable in the form 
P = GnM = pg/2; 

3) Introducing a gravitational mass defect for the point particle in the following way: 
Representing the bare mechanical mass Mq of the point source in the form Mq = 

Jp°° MQ6{r)dr = 47r JJ""" \/—grr{r) p^ij) p{r)dr, one obtains for the mass distribution of 

the point particle the expression p{r) = Mo6{r)/ (47r (^) P^(^) ) = Mo6g{r), where 



5g{r) := 5{r)/ yAny^—grr (r) (r) j is the ID invariant Dirac delta function. The Keple- 
rian gravitational mass M can be calculated using the Tolman formula [5]: 

proc 

M = An p\r)p\r)^i{r)dr = M^^/^J^. (4.6) 



Here we use the relation p' = ^J—gtt 9rr- As a result we reach the relations: ^'^(O) = 

g2vi(o) ^ exp(^-2^^) = < 1 and = GNM\n{^) > 0. (Note that due to 

our convention 0(0) := 1 the component (7ti(r) is a continuous function in the interval 
r G [0, cxd) and gtt{0) = gtt{+0) is a well defined quantity.) 

The ratio Q = ^ = ^/gttW) ^ [O5 1] describes the gravitational mass defect of the 
point particle as a second physical parameter in the problem. The Keplerian mass M and 
the ratio g define completely the solutions ()4.ip . 

Then for the initial constants of the problem one obtains: 



V9i(0) = \ng, ip2{0) = - In- 
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Thus we arrive at the following form of the solutions ()4.ip : 

^2(r) = -lii (yee-'l""" - ee''^"")^ (4.8) 
and the rho-gauge fixing function 

pBRG{r) = Pg{i- Qy^v(y^r/ p(^^ . (4.9) 

An unexpected feature of this two parametric variety of solutions for the gravitational 
field of a point particle is that each solution must be considered only in the domain 
r G [ 0, GnM In (1/^) ) , if we wish to have a monotonic increase of the luminosity variable 
in the interval [po, cxd). 

It is easy to check that away from the source (i.e., for r > 0) the solutions ()4.8|) coincide 
with the HG solution and acquire the well known standard form, when represented using 
the variable p. This means that these solutions strictly respect a generalized Birkhoff 
theorem. Its generalization requires only a justification of the physical domain of variable 
p. In a remarkable accord with Dirac's intuition ^U] the minimal value of the luminosity 
variable for the solutions ()4.8|) is 

po = 2G^M/(1 - g") > pG. (4.10) 

This changes the Gauss theorem and leads to a variety of different important consequences. 
One of them is that we must apply the Birkhoff theorem only in the interval p G [po, 00) ^ 
r G [0,roo). As a result, in this domain all local GR effects like gravitational redshift, 
perihelion shift, deflection of light rays, time-delay of signals, etc., will have their standard 
exact values in gravitational field of the solutions (j4.8j) . 



4.2 Regular Mapping of the interval r G [0, Too) to the whole 
interval r G [0, oo) 

It does not seem to be convenient to work with the unusual radial variable r G [0,roo). 
One can easily overcome this problem using the regular radial gauge transformation 

™ ^ (4.11) 



r/r + 1 

with an arbitrary scale f of the new radial variable r (Note that in the present article 
we are using the same notation r for different radial variables.) This linear fractional 
diffeomorfism does not change the number and the character of the singular points of the 
solutions in the whole compactified complex plane of the variable r. The transforma- 
tion 1)4.111) simply places the point r = Too at infinity: r = oo, at the same time preserving 
the initial place of the origin r = 0. Now the new variable r varies in the standard interval 
r G [0, oo) and the regular solutions ()4.8j) acquire the final form 



(^i(r) = -ln(l/^) 



^ r/r Q / r/r 
r/f+\ \r/f + l^ 

fl ( _^ 1 

V^2(^) ~ ~ ( 2 "^^^^^ ~ Q-^lf^^ 



r/f+1 

(^(r) = 21n(r/r + 1) + ln(r/roo). (4.12) 
The final form of the rho-gauge fixing function reads: 

PPi?G(r) =Pg(i-^^) (4.13) 

The last expression shows that the mathematically admissible interval of values of the 
ratio Q is the open interval (0, 1). This is so, because for = and for = 1 we would have 
impermissible trivial gauge-fixing functions ppijG('^) = l and Pp_rg(?^)=0, respectively. 

The expressions (|4.12p and (|4.13|) still depend on the choice of unit for the new variable 
r. We have to fix the arbitrary scale of this variable in the form r = pc/ ln{l/ g^) = 
G]\fM/ In (^) to ensure validity of the standard asymptotic expansion: gu ^ I — Pg/t" + 
{{pg/^Y) when our last radial variable r goes to infinity. Then the final form of the 4D 
interval, defined by the new regular solutions outside the source (i.e., for r > 0) is: 

ds' = e'^o {de-j^}j -pPRGirf {de' + sin'ed<t>'). (4.14) 

Here we are using a modified (Newton-like) gravitational potential: 

ifJr-M,Mo):= , ^ (4.15) 

' r + GjvM/ln(^) ^ ' 

a coefficient N^{r) = {2(p^)~'^ (e^'^c — 1), and a luminosity variable 

p.„.M^2G„M/(l-.-<=)^:±^l^^if<tl. (4.16) 



These basic formulas describe in a more usual way our regular solutions of Einstein 
equations for r G (0, oo). Therefore we call this radial gauge physical radial gauge (PRG). 

The Eqs. fl4.14|) - ()4.16|) show immediately that in the limit g our solutions tend 
to the Pugachev-Gun'ko-Menzel one Jl], and cp^ (0; M, Mq) = In g ^ —oo. In the limit: 
— > 1 we obtain for any value of the ratio r/pg-. guir/pc, g) — >■ 1, grrirj pc-, q) —1, and 
p{r/pc, g) — > oo. Because of the last result the 4D geometry does not have a meaningful 
limit when g ^ 1. In this case <y9g (r; M, Mq) — > at all 3D space points. Physically this 
means that solutions with point source without mass defect are not admissible in GR. 

5 Total energy of a point source and its gravitational 
field 

In the problem at hand we have an extreme example of an "island universe". In it a 
privileged reference system and a well defined global time exist. It is well known that 
under these conditions the energy of the gravitational field can be defined unambiguously 
[5J. Moreover, we can calculate the total energy of the aggregate of a mechanical particle 
and its gravitational field in a canonical way. Indeed, the canonical procedure produces a 
total Hamilton density "Ktot = ^a=i,2;fi=t,r < V^a,/. - ^tot = 2^ (-pVi^ + pV2^ - e^"^") + 
MQe'^^6{r). Using the constraint (j4.4j) and the first of the relations (j4.7j) . one immediately 
obtains for the total energy of the GR universe with one point particle in it: 

POO 

Etot= / ^totdr = M = gMo < Mo . (5.1) 
Jo 

This result completely agrees with the strong equivalence principle of GR. The energy 
of the gravitational field, created by a point particle is the negative quantity: Egji = 
Etot -Eo = M-Mo = -Mo(l -g)<0. 

The above consideration gives a clear physical explanation of the gravitational mass 
defect for a point particle. 

6 Invariants of the Riemann Tensor 

To reach a coordinate independent description of the geometry of space-time manifold with 
metric (j3.1|) we calculate the invariants of the Riemann curvature tensor. Starting from 
the general consideration in J2] we find in BRG {(f = 0) the following four independent 
invariants: 

I, := e2(^i-2^2)y,;', 

/3 := e2(^-2^^) {-^[^ + ^'^^ - e^^'^/p") /2, 
/4:=e^(--^-^)^;(^;-^'2)/2 

These are linear with respect to the second derivatives of the functions V5i,2 ^ a property, 
which is of critical importance when we have to work with distributions <^i2- Fo'^ the 



(6.1) 



regular solutions ()4.8j) one obtains: 



^Pg Q \PgJ 2 



e{r/pG)-l (._2 4r/pa)^-PG 



2 e(r/p, 



/4 = ^^ (l-,V^/-)^ = ^^^. (6.2) 

The invariants /i,...,4 of the Riemann tensor are well defined distributions. This con- 
firms the general expectations, described in the articles ^3], where one can find a correct 
mathematical treatment of distribution-valued curvature tensors in GR. 

As we see, the manifold Wl^^'^^{g^^{x)} for our regular solutions has a definite geomet- 
rical singularity at the point, where the physical source - the massive point particle is 
placed. The fractional-linear transformation ()4.11|) does not change the character of this 
singularity and does not add new ones. 



7 Global Analytical Properties of the Solutions in 
the Complex Plane of the Radial Variable 

To simplify the notations we introduce dimensionless variables: ( = r/Gj^M - instead of 
the radial variable r in BRG, z = (t/GnM) ln(Mo/M) - instead of the radial variable r 
in PRC, and w = p/ pc - instead of the luminosity variable p. In this Section we shall 
consider these variables as complex valued ones in the corresponding complex plains C,^, 
Cz, C^, and study the basic analytical properties of the corresponding functions. Then, 
up to inessential constant factors, the independent curvature invariants (j6.2p read: 

/i~5(C), /3~(e(C)-l)(l-f^V^)\ /4~e(C)(l-f^V^)' -inBRCor (7.1) 



2\^, 2 



I,^6{z), h^[e\^-^j-lj , I^^e\^-—j l^l-g-j -in PRC. (7.2) 

As we see, the BRG-invariants have a unique analytical property: their only singularity 
in is the center C: C = 0- the compactified complex plain these invariants have 
an additional essentially singular point C = oo, where the jump of the luminosity variable 
w{() = (l — f)^e^'>) on the real axes equals 1. As seen from Eq. ()7.2p . the transition to 
PRG, defined by the fractional linear mapping ()4.1H) . which now reads 

a^) = {\nl/g)^ ^ z{0= f (7.3) 
2+1 [\nl/g) - C 

simply translates this essentially singular point at the position z = —1. From analytical 
point of view the existence of such essentially singular point is the basic difference between 
the descriptions of massive point source in GR and in Newton gravity. 

It is obvious that transition to any other radial variable z, defined by nonlinear trans- 
formation z = z{(), more general than a fractional linear one: z{() = will create 



new singular points in the plane C^, or will change the character of the existing ones, thus 
changing the mathematical and the physical properties of the solution in the plane. 

Hence, from analytical point of view the PRG is a unique radial gauge which defines a 
preferable radial variable r in the problem at hand, satisfying simultaneously the following 
two conditions: 

i) In the physical domain it varies in the natural interval r G [0, oo) and has no other 
singularities in this interval than the very point source. 

ii) The only singular points of the solution in the compactified complex plain are: 
the place of the source at r = and the unavoidable in GR essentially singular point at 
r = — 1, which is placed in the nonphysical domain. 

This important result solves the longstanding problem of the choice of radial variable 
r for point source in GR on a clear theoretical basis. 

The only independent curvature invariant for the HG form of Schwarzschild solution 

()2.5|) is {Ia)hg ~ '^'^ ~ {^ — Q^^ The comparison with Eqs. ()7.2p makes it clear 
that from geometrical point of view the HG solution of EE ()2.5|) is a vacuum one and 
essentially differs from the massive point particle solutions. 

Now we are ready to describe the singular character of the coordinate transition from 
the Hilbert form of Schwarzschild solution to the regular one ()4.12|) . considered outside the 
source. Eq. ()4.13p shows that in this domain the change of the coordinates is described, 
in both directions, by the functions: 

w{z) = (i-qM~' ^ z{w)= ^^^yf^ -1, ge{OA). (7.4) 
V / Inw — m[w — I) 

The function w{z) is regular at the place of the point source z = 0; it has a simple 
pole at 2; = oo and an essentially singular point at z = —1. At the same time the inverse 
function z{w) has a logarithmic branch points both at the HG center of symmetry w = 
and the event horizon w = 1. Thus we see how one produces the HG singularities at 
p = and at p = p^, starting from a regular solution. The derivative 



dz/dw 



ln(l/^2 



w{w — 1) {Inw — ln{w — 1))^ 



approaches infinity at these two points, hence the singular character of the change of 
the variables in the whole complex domain. The restriction of the change of the radial 
variables on the corresponding physical interval outside the source: z G (0, 00) ^ w & 
(1/(1 — g^), 00), is a regular one. 



8 Conclusion 

In the present article we have studied a new, two parameter class of solutions of Einstein 
equations. These static spherically symmetric solutions describe the gravitational field of 
massive point particle with bare mass Mq > and Keplerian mass M (0 < M < Mq). The 
difference between these masses, or their ratio g = M/Mq G (0, 1), define the gravitational 
mass defect for the point particle. Such mass defect was not considered and studied until 
now, because for the standard Hilbert form ()2.5|1 of the Schwarzschild solution "the bare 
rest-mass density is never even introduced" jjlj correctly. 



The new solutions form a two parameter family of metrics on singular manifolds 
M^^'^^{g^i,}, described in details in the present article. 

We have shown the principal role of the massive point source of gravity. It presents 
a natural cutting factor for the physical values of the luminosity variable p G [po,oo), 
where po > pc ()4.1()j) . This happens because the infinite mass density of the matter point 
changes drastically the geometry of the space-time around it. 

Similar geometry of space-time with po = Pg was discovered at first in the original 
article by Schwarzschild j2]. According to Eq. ()4.10|) . such limiting value of the luminosity 
variable corresponds to zero value q = Q oi mass defect ratio. For finite value of M 
this is possible only if Mq = oo. In this sense our work is a proper extension of the 
Schwarzschild one to the physically and mathematically admissible values of the mass 
defect ratio g G (0, 1). 

In full accord with Dirac's suggestion |10^ our cutting of the domain of luminosity 
variable places the event horizon in the nonphysical domain of the variables. This effect 
is well known from the solutions of Einstein equations with massive matter sources of 
finite dimension. 

The mathematical and the physical properties of the new solutions are essentially 
different in comparison with the well known other spherically symmetric static solutions 
to the Einstein equations. All they have different type of singularities at the center of the 
symmetry, which is surrounded by empty space. The previously known solutions were 
often erroneously considered as a solutions for describing of single point mass in GR. 

It is clear that our solutions in generalized functions define in mathematical sense 
the fundamental solutions of Einstein equations, which are complete analogous to the 
fundamental solutions of Poisson equation in Newton theory of gravity. Thus the problem, 
formulated by Feynman in |T| is solved. 

Further study of the new solutions and their physical applications will be given in 
separate articles. 
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